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Abstract The purpose of this paper is to introduce a new way to inquire the quan¬ 
tum cosmology for a certain gravitational theory. Normally, the quantum cosmological 
model is introduced as the minisuperspace theory which is obtained by reducing the su¬ 
perspace where the Wheel-DeWitt equation is defined on using the symmetry provided 
by cosmological principle. Unlike that, the key of our approach is to reinterpret the 
cosmology in a classical dynamical way using a point-like Lagrangian and then quan¬ 
tize the point-like model. We apply the method into Einstein gravity, gravity with 
a cosmological constant and the /(R)-gravity, and get their wave equations respec¬ 
tively. By analsysing the exact solution for the quantum cosmology with and without 
a cosmological constant we demonstrate that the cosmological constant is essential 
and being a tiny positive number. We also show the possibility of explaining inflation 
under the quantum version of cosmology. 


PACS. 04.20.Fy 04.50.Kd 04.60.Ds 

1 Introduction 

The motive of studying a quantum cosmological theory mainly emerges from two parts. First of all, 
classical gravity theory fails in precisely describing topics like the very early universe. To understand 
those tiny scale objects with huge energy, a successful quantum gravity theory is requisite and then 
we can apply it into the cosmological case. However since no such theory is available till now, 
one considerable effort is to view the quantum cosmology as an effective theory which can be 
approximately obtained by modifying the classical theory. 

In classical field theory, the Lagrangian density is a function of one or more fields and their 
derivatives. Then one gets the action 



( 1 ) 


According to the principle of least action, one can get corresponding equations like Maxwell equa¬ 
tions, Klein-Gordon equation, Dirac equation etc. Similarly, in the general relativity, the Einstein- 
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Hillbert actions discribed as 


Seh 


Cd^x 


^—gRdt^x 


(2) 


which can get the Einstein field equation. This Lagrangian formulation leads to get Hamiltionian 
formulation, through Legendre transformation, which would make the transition to quantum me¬ 
chanics easier. An attempt following this idea is called the minisuperspace quantization which is 
first purposed by DeWitt [T]. At that time Wheeler had introduced the idea of superspace, the 
space of all three-geometries as the arena in which geometrodynamics develops, a particular four- 
geometries being represented by a trajectory in this space. Misner had just finished applying the 
Hamiltonian formulation of gravity, developed in the late 1950s and early 1960s, to cosmological 
models with an eye towards quantization of these cosmologies as model theories of general relativ¬ 
ity m, which is another motive for studying quantum cosmology. He invented “minisuperspace” 
and “minisuperspace quantization” or “quantum cosmology”, to describe the evolution of cosmo¬ 
logical spacetimes as trajectories in the finite-dimensional sector of superspace related to the finite 
number of parameters, which describe t =const slices of the models, and the quantum version of 
such models, respectively. 

For the classical cosmology, its dynamic equation is the Friedman equation which is derived 
from Einstein field equation. Under standard model, for the quantum cosmology, a quantum gravity 
theory might be an effective approach. Cosmological minisuperspaces and their quantum versions 
were extensively studied in the early 1970s, but interest in them waned after about 1975 and little 
new work was done until Hawking revived the field in the 1980s [2113], emphasizing path-integral 
approaches. This started a lively resurgence of interest in minisuperspace quantization till now. For 
a quick review of the minisuperspace quantization see Halliwell 0017]. 

In this paper, we would like to introduce a new technique in acquiring the quantum system 
of cosmology via the point-like Lagrangian under a certain gravitational model. According the 
cosmological principle, if we use the spacial symmetry of our univese first, a point-like Langrangian 
can be got as bellow 

L = J CdQ = J ./^Rdn , (3) 

where g is the determinant of Friedmann-Robertson-Walker(FRW) metric, R is the scalar curvature 
and dfi represents the spacial volume element. The Ensitein-Hillbert action can be written as 




/ 


L dt 


(4) 


here t is the cosmological time, which can be treated as a classical mechanics action. Then quantize 
this model through the classical canonical quantization. It is much easier than the normal process 
but reaches the similar result. We will then apply this technique into a more complicated case, 
which contains a cosmological constant, and even more general, for the f{R) gravity. We will give 
the exact solution to the quantum universe with a cosmological constant and discuss its meaning. 


2 The Friedmann Equations 


In the normal sense, we get the cosmological model under a certain kind of gravity by applying the 
cosmological principle into the gravitational field equation which is obtained from variation of the 
action. That is to say, considering a action with matter field 


S= / ^Rd^ 


■ Sm 


(5) 
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by varying the action eq.® with respect to one can get the Enistein Field Equation. After 
setting the metric in the equation to the FRW metric 


ds^ 


—dt d (t) 


dr^ 


1 — 


+ {d9^ +sm‘^ edip^) 


( 6 ) 


and the energy-momentum tensor with a perfect fluid 

= {P + pW^U'' + 


(7) 


it would end up with the Friedman Equations. 

Surprisingly, one can apply the cosmology principle directly into the action before variation. 
This modifies the Lagrangian to a classical point-like one with respect to the scale factor. By the 
variation of the point-like action, it will give exactly the same equations as those in the normal 
process. Such method has been extensively used in studying scalar field cosmology [snanQ] , non- 
minimally coupled cosmology scalar-tensor theory [13], multiple scalar fields [H], vector 

field [15], Fermion field [TB], /(T) gravity [T7], f{R) gravity [T5], high order gravity theory [T^ . 
GaussBonnet gravity [20] and so on. Without loss of generality, we get the action as below by 
setting K = 0, 

S= —6aa^dt + Sm , (8) 


the variation with respect to a and take note of SSm/^o, = SSM/Sg^i, ■ 5g^^/6a 
one of the fridemann equation 

d 

2- + — = -p 

a a 

By combining the conservation equation of perfect fluid p = —3(p-|-p)d/a, all Friedmann equations 
can be derived from this point-like Lagrangian. 

Moreover, this point-like Lagrangian may lead us to quantiza the cosmology, one would get 
an easier wave equation and its analytic solutions. In addition, through the classical canonical 
quantization, the evolution of the cosmology which can be discribed by the scale factor, may be 
different from the classical cosmology as it has the quantum effects. Next, we will apply this method 
into the quantum cosmology. 


= 6p/a, one can get 
(9) 


3 Point-Like Quantum Cosmology 

For simplicity, let us restrict our discussion to the flat FRW universe(K = 0) in the case of without 
considering the matter field. Recall the point-like Lagrangian for a flat FRW universe under Einstein 
gravity, its effective part is 

L = -6aa^ . (10) 

Its canonical momentum is 

dL 

TTa = -^ = -12ad , (11) 

from which we can get its Hamiltonian 

7r2 

H = TTad — L = —6ad^ = — — 

24a 


(12) 
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Following the process of standard canonical quantization, according to the Weyl rule m , we can 
get the Hamiltonian operator 


H =+ 27ra-Tra + -nl) , (13) 

24 4 a a a 

then we should replace the canonical momentum tt^ by the operator of momentum —ida and get 
the wave equation of the scale factor a 


.dip 11 d'^ip 1 dip 1 
^ dt 24 a da^ a? da 2a^ 


(14) 


Since the scale factor a is a non-negative real number in cosmology, this wavefunction should be a 
function defined only on the right half of the real line with its value being a complex number at 
any given time. 

Assuming ip{a,t) takes the form of e“*'^‘(^(e; a) and pluging it into ea. lfldll . we get the eigen 
equation 

^[^cP"{a) - ^(p'(a) + ^(p{a)] = e(p{a) , a^O . (15) 

The solution of this equation depends on whether the eigenvalue is positive or not. 

When e > 0 we can first rescale the variable ^ = (24e)^/^a to drop the parameters in eq. (USD 
and modify it to 

^^cP"-^cP' + {^-e)cP = 0 , ^^0 . (16) 

Through a variable substitution of z = and introducing a new function u{z) with the 

relationship (p{z) = ^ u{z), we get 


u\z) + -u\z) - (1 + (v^/^) _ Q ^ 

which is a modified Bessel equation whose solutions can be described by a linear combination of 
the hrst modified Bessel function I^/^{z) and the second modified Bessel function K^^^{z). 

However increases in exponential form when z goes to infinity, which shows great 

divergent trend and can not be accepted as a reasonable wavefunction. Therefore, recovering all the 
transformations we made, we get the eigen function with a given eigenvalue described by ea. lfTSl) of 


cp+{e; a) = . (18) 

The asymptotic expansion of the second modified Bessel function with huge arguments can be 
described as 



which indicates our wavefunction ip'^ behaves like 

<))+(e;a) - , (20) 


and will decay to zero at an extremely quick rate with a getting huge enough. This means is 
normalizable and is suitable for being a wavefunction. 
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For tiny arguments, the second modified Bessel function has the following asymptotic form, 

Applying this into the wavefunction eq. m, one can find when a —>■ 0, the wavefunction 

^+(e;a) - (22) 

goes zero which can avoid the possibility of a cosmological singularity. 

When e < 0, in order to keep ^ being positive, we should rescale the variable in form of ^ = 
(—24e)^/^a. Now equation eq. m becomes 

+ = 0 . (23) 

Again we apply the substitution z = 2^^/^/3 and let ^ u, then we end up with a Bessel equation 

u'\z) + -u'{z) + (1 - = 0 . (24) 

z z^ 

Normally we choose the Bessel function of first kind J^j^iz) as one of the basis for its solution 
space. The other base can have different choice among which the most convenient should be the 
Bessel function of second kind also known as the Weber function or the Neumann function. 

So for every distinct eigenvalue, the quantum system has the degenerate degree of two with its 
independence basis chosen as 


/'i(e;a)-( 24e)3aJv2( a^) 

3 0 

(25) 

/'2(e;«) = ( 2Ae)3aY^{ ~ 02 ) 

3 0 

(26) 


For great arguments, the Bessel functions behave like the following 

T 


W ^ \/ £: (cos (^ - Y “ 0 + ^)) 

Therefore, when a goes to infinity, the wavefunctions have asymptotic expansions 


(j>^ (e; a) ^ cos 


f A^/^ 3 2v^-f3 


■ 02 — 


12 


, 1 . 3 2^2 + 3 

02 (e; a) ^ a sin | —— 02 - 1 




12 


(27) 

(28) 


+ 0{a-i) , 

(29) 

+ 0{a~i) 

(30) 


that somehow behaves like sine-cosine with a raise to the power of 1/4. thus 0^ , 02 are divergent 
when z goes to infinity, which cannot be accepted as a reasonable wavefunction. 
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As shown from our discussion, only (jy^ can act as the wavefuntion and is capable of avoiding 
the cosmological singularity. Its eigenvalue e would be a continue spectrum which can has value 
from zero to infinity. Normally, for any initial states C£|e)(here |e) is the nomalized eigen state), 
the expectation value of a? evolves like 

aP-{t) = '^CeCe'{e\a^\e)cos{{e - e)t + 0e,e')+ '^\ce\^{e\a^\() ■ (31) 


Therefore, the evolution of can be any functions with specific initial state, including the inflation 
of the early universe. However, many theory consider the cosmology constant is essential, like 
ACDM model and Quantum Field Theory. So we would like to consider a quantum cosmological 
model under some more generalized gravitational models like gravity with a cosmological constant. 


4 Quantum Cosmology with A Cosmological Constant 


The point-like Lagrangian for flat FRW universe under the gravitational model with a cosmological 
constant has the form 

L = -6aa^ + Aa^ . (32) 

The additional term of AaA does not contain the derivative of a, thus will not change the form of 
canonical momentum we get in eq. (HU. So the Hamiltonian in this case is 


H = - Aa^ 

24a 


(33) 


The same quantization method we can get the wave equation 



1 1 gy 

24 a da^ 


1 dij) 


+ _^)_HaV 


(34) 


Its eigen equation is given by 


U' + ( A - 24Aa4 

a 2a^ 


24ea)(^ = 0 


(35) 


Next we will consider the asymptotic behavior of the equation. 

When a —?> 0, the wave equation can be described as 

(f)" - -cj)' + = 0 , (36) 

a 2a^ 

obviously, 

. (37) 

When a —)> oo, the equation behaves like 

(j)" - {2AAa^ + 2Aea)(j) = 0 , (38) 


it has a special solution with the form like 


(39) 
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Generally, assuming takes the forms like 


4>{a) = a 


_ „l + V2/2 . 


■ u{a) 


(40) 


or 


0(a) = a} ■ u{a) . (41) 

Applying ea. (l40l) into the eigen equation eg. (1351) and setting 9i^‘^ = 24A, one would get the equation 
for u{a) as bellow 

au"{a) -\- ( 6 ^a^ -I- 1 -I- '/2)u{a) + [(9 -I- 3'/2)iy — 24e]a^tt(a) = 0 . (42) 

With substitution z = —2va?, rearrange ea. d^ 

Q I . A) /I, 

(43) 


zu''{z) + - z)u'{z) - = 0 


which is a confluent hypergeometric equation that is also known as Rummer’s equation. It has a 
solution described by Rummer’s function defined as 


F(a,7,z) = 


n—0 


j{n) 


(44) 


where in our case a = (3 -I- •\/2)/6 — 4e/3z^ and 7 = (3 -I- •\/2)/3. Here the symbol x^'^\x = a, 7 ) 
stands for a rising factorial defined as 

= 1 , (45) 

= a;(x-I-1) • • • (x-I-n — 1), n>l . (46) 

Since 7 is not a integer, another solution independent with ea. (H^ can be simply introduced as 


{z) = z^ '*'F(a - 7 -I- 1 , 2 - 7 , z) 


(47) 


Recovering from the substitution we made and considering v can be either positive or negative, 
it seems that we will get four independent eigen functions for any given eigenvalue e 


,3-|-'\/2 4e 3 -|-3 

6 ^3.' 3 ' 

(48) 

A-\/2 4e 3-\/2 „ 

6 ^3.- 3 ' 

(49) 


0^(e;a) = a^ \/2/2g±t 

where v = 2\/ 6A/3. However Rummer’s function obeys the property of Rummer’s transformation 

F(a,j,z) = e^F(-y - a,-f,-z) . (50) 

So in fact one can check that (j)i 2 — 4>i2i there are only two independent solutions for each 
e. In addtion, when apply eg. (HD) into eg. (l55l) we get the same two solutions except a constant 
coefficient according to Rummer transfermation. 

For great arguments. Rummer’s function can be approximately expanded as 


F{a,j,z) - r( 7 ) 




g2^Q;-7 


F{a) F (7 — a) 


(51) 
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Applying it into ea. (H51) and ea. lH^ . we get the asymptotic expansions of these wavefunctions when 
a goes to infinity, 


where 


01 - 

i + V2 3 + 

'F){ g , v) + e{ g ,v) 

(52) 

02 

!?> — \f2 3 — 

g , i/) + 0( g , 1 /) 

(53) 


“ (2iz)(«-4«^/30r(^ + 4e/3i/) 

(54) 


Since v = 2\/KA/‘i, the sign of A will decide whether v is real, and therefore decide how the 
asymptotic expansions behave. 

When A < Q, V = is an imaginary number and makes the modulus of both 

exp(±zza^) and being unit for any real a. So the asymptotic expansions can be simplified to 


(/' 1.2 ^ 0{a 


(55) 


This descending with the order of minus one half is too slow, thus neither of the two eigen functions 
is normalizable. 

When A > 0, with v being a real number, the behavior of O is completely determined by 
the exponential term exp(±zza^). Fortunately, we know that r{z) explodes at non-positive integer 
points, thus a carefully selected eigenvalue e can make the exploded term 0 (^, v) in the expansions 
vanish. 

For (j)i it requires 


(6n -l- 3 + \/2) rA 


el.' = 


n = 0 , 1 ,... 


(56) 


And we have the fact that for F(a, 7 + 1 , z) whose a is a non-positive integer, it can be described 
by Laguerre function 

L^uHz) ■■= . (57) 

So the eigen states for eigenvalue after normalization is written as 


where Cn'^ is the normalization factor. 

Same discussion also holds for solution 02 , giving its eigenvalue 


e(2) = - 


( 6 n -I- 3 — ^/2) l~A 

2 V "6 


u = 0,1,... 


and the normalized eigen state 


(58) 


(59) 


(60) 


( 2 ) 

where Cn ^ is its normalization factor. 
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From the definition of Kummer’s function eg. (1441) we can see that for tiny arguments, 


^i(e ; a) ^ — i/a^)(l-I-2z/—a^) , (61) 

7 

02 (e ;a) ^ . (62) 

2-7 


Therefore when a goes to zero, (pi and p 2 are all vanish, hence the cosmological singularity is 
naturally avoided. Thus we can choose the two sets of eigen states to be the basis of the quantum 
system: en \ and en\ 

For any real number a, the first two Laguerre polynomials are 


4^) = 1 , (63) 

= i + j-x . (64) 


So the first set of eigen states with eigenvalue of the highest two are 




(65) 

( 66 ) 


If a certain state is a combination of only these two states \ip) = -I- at initial (|ap + 

= 1), it will evolve with respect to the cosmological time t as 

= + • (67) 

The evolution of the average measurement can be calculated out: 


a2(t) ={il;,t\a^\pj,t) 

=|a|2(0|a2|0) + a*/3(0|a2|l)e*(^o'’-^™)‘ + 
=25i(a*/3e*(^o'’-^i'’>‘)(0|a2|l) + (|a|2(0|a2|0) + |/3|2(l|a2|l)) 


=2co,i|a||/3| cos 


'Vm 


t-e] + (co.olap + ci,i|/3p) 


where Cn,m = 9 = Arg(a) — Arg(/3). Fortunately, a certain state which is a combina¬ 
tion of can get the similar result as {e^n ^ — em^) is also always some integer times of \fKAl2 

according to eq. dSi). 

Therefore it is reasonable to assuming that our universe is in a state of combination of |n^^^) or 
This solution suggests a pulsing universe with a characteristic time of Att^/TJqA. Actually 
although we do not know which state the universe is at a certain time, we can prove that it always 
rebounds with the same characteristic time no matter how the initial coefficients for the state of 
each energy level are given. It is clear that the evolution of a? is always described by 

Q?'(ip = ^ ( 'l]nm COs((£:n Sm)i 4“ 9nm') 

n,m ^0 


(68) 
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We can see £„ — £m is always some integer times of which ensures that a(t) has a period of 

47ry^l/6yl. 

Considering the universe as we observed is experiencing an accelerating expansion now, it is 
reasonable to assume it is still in the first quarter of the period, implying the cosmological constant 
A should be no bigger than tt^/GTq where Tq stands for the cosmological time till now. 

Another interesting thing is that, observing the combination of solutions can provide a square 
wave, this quantum system must have some special states that may let rise as fast as possible 
at some certain time to- For the simplest case, considering a state composed by eigen states with 
real coefficients lijjA) = Sn 'rr!,e“*^"‘|n), if its scale factor evolves like 

+ C , (69) 

then choose only the coefficients of the first 2N — 1 states to be non-zero and real, we know they 
will satisfy the polynomial system 

A / i_ ( _iy\ 

TnTn+lCn,n+l = J f - J -j , 1 < 1 ^ 2Af - 1 , (70) 

There are in total 2N coefficients needed to be fixed. The normalization condition together 
with eq. dZOl) exactly give the same number of equations from which the coefficients can be solved. 
Therefore we can satisfy ea. (p511 for any N as large as we wish. That gives the possibility of an 
expansion of a at any velocity, which may generate an inflation with the speed even faster than 
exponential level as normal understanding. 


N 




COS 


'Vm 


(2fc- 1) 


A;=0 


5 Quantum Cosmology of f{R) Gravity 


For more general cases, we consider the quantum model of a flat FRW universe under /(i?) gravity. 
More detailed discussion on the point-like model of f{R) universe can be found in the works of 
Capozziello [18] . The point-like Lagrangian with no term of matter will be like 


L = {f - fRR)a^ - QfRRRa?a - 6 /^ 00 ^ 

(71) 

and the canonical momentums for a and R are respectively 


TTa = -QfRRRa^ - l2fRaa , 

(72) 

T^R = -GfRRO^d 

(73) 

So the canonical energy is 


El = -{f - fRR)a^ - efRRRa^a^ - QfRaa^ 

(74) 

From eg. (1751) we can directly read that 


QfRRa^ 

(75) 


Plug it into ea. (|7^ and get 

QfRRRa'^ = -ttq + --f^TTR 
a JRR 


(76) 
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Applying them to ea. dTH) . the Hamiltonian for the system becomes 


H =-if - fnR)a^ + 
= -if-fRR)a^- 


^IrroP' 
1 1 
6a2 fRu 


T^RT^a 


2 fR 
a fRR 

1 Ir 

6a3 P 


TTR 


RR 


So the wave equation that describes this quantum system is 



\ 2 

^Irro^ J 


■9 r 


1 Ir 52 
6a3/|^ dR^ 




Irrr d _ 
PRR da 



1 1 


52 




6a2 fnn dadR 
Irrr ^/rrr 


mfi 


ff 


RR 


P 


RR 


1 fQfRfRRR 


■I RR 


24a3 V 


7 

Jrr 


OR 




P 


RR 




\f - IrR) 


If' 


1 

12 a2 


a, R) is a function of cosmological time t, scale factor a and Ricci scalar R. 

The equation relies on the form of f{R) to be exactly solved. However as a linear partial differ¬ 
ential equation, its coefficients of all the second order terms satisfy the fact that A = l/(144a'^/|.^) 
is positive on the whole a — R plane. Therefore the eigen equation of the operator H is a hyperbolic 
equation and can be transformed into a wave equation. 

We need to point out that a and R have been separated via Palatini formalism. In this case, 
their relation is linked by one of the equations of motion rather than a given definition. So after 
the quantization, this relation must have degenerated to be statistically satisfied. That means even 
a flat universe of small scale or a huge scale universe with large curvature which are not normally 
allowed in the classical case will also have contribution to the possibility. 


6 Conclusion And Discussions 

The purpose of this paper is to introduce a new approach to inquire the minisuperspace model 
without seeking the Wheeler-DeWitt equation for a certain gravitational theory. 

The technique is to apply the cosmological principle directly to the action of a gravitational 
system before variation, and reform the Lagrangian of geometry to a classical point-like one. It is 
obvious that such a process of taking the metric of a cosmological model which is truncated by an 
enormous degree of imposed symmetry and simply plugging it into a quantization procedure should 
not give an answer that is in any way an exact solution. However, strange enough, we have seen 
that the variation of this point-like Lagrangian gives the right equation of motion (the Friedmann 
equation) to describe the universe. 

By quantizing this semi-classical system described by the point-like Lagrangian, we represent 
a quantum system that is very similar to the minisuperspace from reducing the superspace where 
the Wheel-DeWitt equation is defined on. The only difference is, in our situation, for solving a 
semi-classical Schrodinger equation we need the concept of the eigenvalue e of the Hamiltonian of 
the system which does not exist in the classical minisuperspace theory. 

It is very natural to apply our technique beyond the Einstein gravity to the gravitational model 
with cosmological constant and more general f{R) gravity with the help of Palatini formalism and 
respectively get their quantum cosmological model. This especially opens the gate for considering 
quantum cosmology of f{R) gravity. 

As the second aspect of our work shown in this paper, we give the exact solutions of the quantum 
systems we get under Einstein gravity with and without cosmological constant. We find that the 
existence of a tiny positive cosmological constant is reasonable 
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We prove that all possible states in such a legal quantum cosmological model predict pulsing 
universe with the same period of a cosmological characteristic time that is inversely proportional to 
the square root of the cosmological constant. Considering the enormous amount of time the universe 
has existed, the cosmological constant must be extremely tiny. 

Moreover, we show that this quantum system contains states that allow expansion at any speed 
as fast as possible, which could probably provide a motivation for inflation. 
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